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The universe might have had a prolonged exponentially expanding phase caused by its being stuck in a metastable state 
of the grand unified phase transition. The only way that it could exit from this exponential expansion without introducing 
too much inhomogeneity or spatial curvature would be through a homogeneous "bubble" solution in which quantum 
tunnelling occured everywhere at the same time. This would produce more baryons than the conventional scenarios. 



There has been considerable interest recently in 
the possibility that the universe might have under- 
gone an exponentially expanding "inflationary" stage 
caused by the vacuum energy of a supercooled meta- 
stable state of a first-order phase transition associated 
with the breakdown of the grand unification symme- 
try [1-4] . Such an exponential expansion might pro- 
vide the answer to a number of problems in cosmolo- 
gy such as 

(1) Why are there so few magnetic monopoles 
[1,5]. 

(2) Why is the universe so spatially homogeneous 
and isotropic? 

(3) Why is the universe so nearly flat, i.e. why is it 
expanding at almost exactly the critical rate to avoid 
recollapse? 

The main difficulty with this inflationary scenario 
has been to find some way of getting out of the ex- 
ponentially expanding stage without creating large- 
scale inhomogeneities that would be incompatible 
with observation. In his original paper Guth [2] sug- 
gested that small bubbles of the broken symmetry 
phase would form by vacuum tunelling [6] and then 
expand by a large factor before colliding with other 
bubbles. However, as he himself admitted, and others 
have confirmed [7] , one would end up with a very 
inhomogeneous universe dominated by a few very 
large bubbles. 

Most of the papers that have discussed the exit 
from the inflationary stage have treated it essen- 



tially as a phase transition problem in flat spacetime 
and have neglected the curvature and finite horizon 
size of the universe. The aim of this letter is to show 
that these effects play a very important role and can 
lead to a phase transition which occurs simultaneous- 
ly at all points of space thus not creating any inhomo- 
geneity or magnetic monopoles. If the adjustable 
parameters in the effective potential for the Higgs 
fields lie in a certain range, the spatial flatness of 
the universe would also be accounted for. 

The effective potential of the Higgs fields has 
the general form 

V(<t>) = ?(m 2 +£R +cT 2 yp 2 

+ ^a 2 0 4 (log0 2 /0^ - \ +m 2 /a 2 <pl) 

+ | a 2 0 4_i w 2 0 2 f (1) 

where 0 is some measure of the magnitude of the 
Higgs field and 0 O is its expectation value at zero 
temperature in flat space with V(<p Q ) = 0. The quan- 
ties m 2 and % are renormalised parameters which 
can be given any value. % = \ is a natural choice be- 
cause it corresponds to conformal invariance when 
m = 0 and it is not renormalized at one loop though 
it is at two loops. The quantity c is a numerical fac- 
tor of order 1 which depends on the group and the 
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representation. Neglecting small corrections arising 
from the %R term, the energy density of the universe 
will be 

p = (77 2 /3OMr>7 4 + H0^), (2) 

where g(T) is the effective number of spin states of 
relativistic particles. The rate of expansion will be 
given by 

H 2 = (87r/3m 2 )p - k/S 2 (3) 

where H = 5/5 and 5 is the scale factor of the uni- 
verse. 

At very high temperatures the universe will be in 
the symmetric phase 0 = 0. If m 2 +%R> 0, the sym- 
metric phase will remain a local minimum of the 
potential as the temperature falls although it will 
cease to be a global minimum at some critical tempe- 
rature T c of the order of a0 Q . However, in realistic 
models the probability to tunnel to the broken sym- 
metry phase with 0 ~ 0rj is very low [8] . The uni- 
verse will therefore continue to expand in the sym- 
metric phase until the thermal energy, the first 
term in eq. (2), becomes small compared to the 
vacuum energy V(0, T). This will lead to an ex- 
ponential expansion and the universe will rapidly ap- 
proach a de Sitter state which depends only on the 
value of V(0) through the relation 

// = [8i:V(0)l3m 2 ] V 2 (4) 

but which is otherwise independent of the initial 
conditions and the initial value of k. This is very 
similar to the way that a gravitational collapse rapid- 
ly approaches a stationary black-hole state which 
depends only on the mass and angular momentum 
but which is otherwise independent of the nature 
of the collapsing body. Following the black-hole 
analogy, we shall therefore assume a cosmological 
"no hair" theorem and treat the universe as a de 
Sitter model. 

In the gravitational collapse case , the black hole 
has an effective temperature [9] T B = m 2 (8irM)~ l . 
Similarly, de Sitter space has an effective temperature 
[10] T s = H/2ir. In the case of grand unified phase 
transitions this will be of the order of 10 11 GeV. In a 
de Sitter universe it is impossible to define a tempera- 
ture smaller than T s . It is therefore incorrect to use in 
the effective potential (1) the redshifted temperature 
of the universe when this falls below T„. In fact one 



18 March 1982 

can regard the %R term as representing the contribution 
of T s to the effective potential. 

In accordance with our viewpoint that de Sitter 
space has no hair, we shall calculate the probability of 
tunnelling from the local minimum at 0 = 0 to the glob- 
al minimum by looking for solutions of the classical 
coupled scalar and gravitational field equations that 
are near the euclidean version of de Sitter space, i.e. 
a four-sphere of radius H r 1 . 

If m 2 > 2// 2 (l-6£), there are inhomogeneous 
"bubble" solutions as there are in the decay of the 
false vacuum in flat spacetime [6] . Gravity and the 
curvature of the universe do not have such effect ex- 
cept to make the bubble slightly smaller [11]. How- 
ever, as mentioned earlier, this case would lead to un- 
acceptably large inhomogeneities in the universe to- 
day. 

If2// 2 (1 -6£)>m 2 >- 12# 2 £,the potential (1) 
would lead to inhomogenous bubble solutions in flat 
spacetime whose radius was greater than that of the 
de Sitter space,// -1 . In this case the only euclidean 
solution of the coupled equations, apart from 0 = 0, 
is the homogeneous solution 0 = 0j where F(0,) is the 
local maximum of V on a four-sphere of radius 
where// 2 = 8irF(0i)/3wjp. The tunnelling probability 
per unit four-volume is of the order of (m 2 + %R) 2 e~ B 
where 

fl=i«}[i/no)-i/n*i)] (5) 

is the difference between the combined gravitational 
and scalar field actions of the 0 = 0j , and the 0 = 0 
solutions. Even though B may be very large and the 
tunnelling probability very small, this does not matter 
because the universe will continue in the essentially 
stationary de Sitter state until it makes a quantum 
transition everywhere to the 0 = 0] solution. 

The 0 = 0j solution is unstable as is reflected by the 
fact that it has a negative perturbative mode 0 = 0] + e 
where e is a constant on the four-sphere. One would 
therefore expect the solution in lorentzian spacetime 
to evolve according to the classical field equations with 
0 running downhill from the maximum at 0 = 0j to 
the global minimum at 0 = 0 Q . The solution will pre- 
sumably be spatially homogeneous but evolving in 
time so it can be represented by a Friedmann— 
Robertson— Walker metric with 0 a functionof time 
only. One can introduce Robertson— Walker coordi- 
nates with K = +1 ,0 or — 1 into de Sitter space though 
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Fig. 1. The scale factor of the universe starting from a k = 
+ 1 time symmetric state with <t> = $i + e. 

the coordinates cover the whole space only in the first 
case. It would seem reasonable to start the classical 
evolution of the solution on a spacelike surface whose 
curvature (positive or negative) is of the order of the 
radius of the de Sitter space. The most natural 
choice would seem to be the surface of time symmetry 
in the K = +1 coordinates but the conclusions in the 
K = —\ case are very similar. 

While the scalar field remains near 0 = 0j , the ex- 
pansion of the universe will continue to be nearly ex- 




Fig. 2. The value of the scalar field in the same universe. 



ponential or, in the K = +1 case, S = c °sh Hyt 
where S is the scale factor. The timescale for 0 to run 
downhill from the maximum at 0 = 0j is 

3^/(m 2 + 12?//?) . (6) 

When 0 approaches the global minimum 0 O it will 
start to perform damped harmonic osciallations around 
0 = 0 O and the expansion of the universe will change 
to that of a model filled with non relativistic matter. 
Figs. 1 and 2 show the numerical solutions for S and 
0in the caseA: = +l,w 2 + \2\h\ = 7AH\. 

One can think of the oscillations of 0 about 0 O as 
corresponding to a very large number of Higgs particles 
in a coherent state. Presumably they will decay into 
light gauge particles and fermions which will have an 
approximately thermal distribution at a temperature 
T 2 given by (rr 2 (30)g(T 2 )T^ * F(0j). The decay of 
the Higgs particles would probably give a much larger 
baryon asymmetry than in the normal scenario. 

The radius of the curvature of the space sections 
of the universe at the time of thermalization will be 
of the order Z/^ 1 cosh //j f 2 where ? 2 is the time from 
the 0 = 0j state to thermalization. If one assumes that 
this is a few times the timescale (6) for 0 to run down- 
hill, one finds that the universe would have been suf- 
ficiently spatially flat to expend to its present size 
with a nearly critical density provided that rn 2 - + \2%H\ 

This seems the only situation in which one can have 
a significant exponential expansion at the GUT era 
without introducing either too much inhomogeneity 
or too much spatial curvature. 

Finally, if m 2 + 12|//j < 0, the symmetric state 
0 = 0 will cease to be a local minimum of the poten- 
tial before the temperature drops to T s so there will 
not be a prolonged exponential expansion phase. 

We are very grateful to A.D. Linde for discussions 
which gave rise to this paper. 
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